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Abstract 

The  evolution  of  microtwins  during  high  temperature  creep  deformation  in  a  y'  strengthened  Ni-based  superalloy  has  been  investi¬ 
gated  through  a  combination  of  creep  testing,  transmission  electron  microscopy  (TEM),  theoretical  modeling,  and  computer  simulation. 
Experimentally,  microtwin  nucleation  sources  were  identified  and  their  evolution  was  tracked  by  characterizing  the  deformation  sub¬ 
structure  at  different  stages  of  creep  deformation.  Deformation  is  highly  localized  around  stress  concentrators  such  as  carbides,  borides 
and  serrated  grain  boundaries,  which  act  as  sources  of  a/2(\  10)  matrix-type  dislocations.  Due  to  fine  channels  between  the  y'  particles, 
coupled  with  a  low  y  matrix  stacking  fault  energy,  the  a/2{l  10)  matrix  dislocations  dissociate  into  a/6(l  12)  Shockley  partials,  which 
were  commonly  observed  to  be  decorrelated  from  one  another,  creating  extended  intrinsic  stacking  faults  in  the  y  matrix.  Microtwins 
are  common  and  form  via  Shockley  partial  dislocations,  cooperatively  shearing  both  the  y  and  y'  phases  on  adjacent  [111}  glide  planes. 
The  TEM  observations  lead  directly  to  an  analysis  of  dislocation-precipitate  interactions.  The  important  processes  of  dislocation  disso¬ 
ciation  and  decorrelation  were  modeled  in  detail  through  phase  field  simulations  and  theoretical  analyses  based  on  Orowan  looping,  pro¬ 
viding  a  comprehensive  insight  into  the  microstructural  features  and  applied  stress  conditions  that  favor  the  microtwinning  deformation 
mode  in  y'  strengthened  Ni-based  superalloys. 
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1.  Introduction 

Ni-based  superalloys  are  a  class  of  alloys  that  are  pre¬ 
dominantly  used  in  the  aerospace  industry  in  hot  sections 
of  gas  turbine  engine  components.  These  alloys  possess 
an  attractive  combination  of  attributes,  such  as  high  tem¬ 
perature  mechanical  properties  (strength,  creep,  and  fati¬ 
gue)  and  resistance  to  high  temperature  oxidation.  The 
high  temperature  mechanical  properties  are  largely  dictated 
by  the  presence  of  ordered,  Ll2  structured,  Ni3Al-based  y' 
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precipitates,  which  are  coherently  embedded  in  a  solid  solu¬ 
tion  y  matrix  (fee  structure).  Since  these  materials  are  uti¬ 
lized  at  elevated  temperatures,  time-dependent  creep 
deformation  is  of  principal  concern.  There  have  been 
numerous  detailed  transmission  electron  microscopy 
(TEM)  studies  focused  on  elucidating  the  creep  deforma¬ 
tion  mechanisms  and  determining  the  associated  implica¬ 
tions  for  creep  rate  controlling  dislocation  processes 
[1-8].  Mechanisms  previously  identified  include  dislocation 
climb,  looping,  and  stacking  fault  related  y'  shearing.  A 
particular  mechanism  that  has  been  reported  more  recently 
following  high  temperature  creep  deformation  is  microt¬ 
winning  [9-17],  which  is  surprising  since  twinning  in  fee 
crystals  is  generally  considered  a  low  temperature,  high 
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strain  rate  deformation  mode,  in  contrast  to  the  high  tem¬ 
perature,  low  strain  rate  creep  conditions  that  have  been 
observed  in  these  alloys  [18]. 

As  with  many  dislocation-based  deformation  mecha¬ 
nisms  commonly  observed  in  y'  strengthened  Ni-based 
superalloys,  dislocation-precipitate  interactions  are  highly 
dependent  upon  alloy  chemistry,  y'  precipitate  size  and 
morphology,  y  channel  width,  crystal  orientation,  stress 
and  temperature.  The  crystallographic  orientation  stress 
dependence  of  microtwinning  in  a  single  crystal  SRR99 
superalloy  was  shown  by  Ardakani  et  al.  [10],  where  it 
was  determined  that  activation  of  { 1 1 1 }  ( 1  1 2)  and 
{ 1 1 1 }  ( 1 1 0)  slip  systems  occurs  when  the  tensile  stress  axis 
is  aligned  with,  or  rotated  towards,  the  [Oil]  and  [111] 
crystallographic  directions.  Similarly,  uniaxial  creep  exper¬ 
iments  of  single  crystal  CMSX-4  by  Knowles  et  al.  [11,12] 
reported  in  greater  detail  how  the  activation  of 
{ 1 1 1 }  ( 1 1 2)  slip  systems  can  initiate  both  stacking  fault 
shear  and  twinning.  It  was  suggested  that  microtwins  form 
by  the  motion  of  at 3(1 12) -type  superlattice  Shockley  par¬ 
tial  on  adjacent  [111]  glide  planes  since  the  Ll2  structure 
of  the  y'  precipitates  would  be  preserved,  as  shown  in 
Fig.  la.  However,  this  would  effectively  require  high  stres¬ 
ses  for  such  a  large  displacement. 

According  to  Christian  and  Mahajan  [19]  deformation- 
induced  twinning  in  Ll2  superlattice  structures  would 
result  in  the  production  of  a  pseud  otwin  characterized  by 
an  imperfect  superlattice  with  incorrect  atomic  ordering  if 
sheared  by  <r//6  ( 1 1 2)  Shockley  partial  dislocations,  as  illus¬ 
trated  in  Fig.  lb.  Atomic  rearrangement  would  thus  be 
required  to  convert  the  pseudotwin  structure  to  a  true  twin 
structure  with  the  correct  nearest  neighbor  atomic  bonds. 
Accordingly,  Kolbe  proposed  a  hypothesis  whereby 


microtwins  in  y'  strengthened  superalloys  form  by  the  pair¬ 
wise  passage  of  identical  a/ 6(1 1 2)  Shockley  partial  disloca¬ 
tions  that  shear  through  both  y  and  y'  on  adjacent  [111] 
glide  planes  [20].  Considering  the  fact  that  shearing  of  the 
Y  precipitates  by  a  single  Shockley  partial  dislocation 
would  form  a  high  energy,  complex  stacking  fault  (CSF), 
and  passage  of  a  second  a/ 6(1 12)  Shockley  partial  would 
create  a  two  layer  CSF,  Kolbe  proposed  that  diffusion- 
mediated  atomic  reordering  in  the  wake  of  the  leading  par¬ 
tial  could  convert  the  two  layer  CSF  into  a  superlattice 
extrinsic  stacking  fault  (SESF)  by  local  rearrangement  of 
atoms  to  eliminate  the  high  energy  Al-Al  first  nearest 
neighbor  bonds.  Although  direct  TEM  evidence  was  not 
provided  at  the  time,  it  was  later  substantiated  through 
direct  TEM  observations  that  microtwins  actually  formed 
when  both  the  y  matrix  and  y'  precipitates  were  sheared 
conservatively  by  identical,  paired  «/ 6(1 12)  Shockley  par¬ 
tial  dislocations  traveling  on  adjacent  {111}  type  octahe¬ 
dral  glide  planes  [14,15]. 

In  a  more  recent  study  the  concept  of  reorder-mediated 
shearing  of  y'  precipitates  was  conclusively  identified  as 
playing  a  role  in  microtwin  formation  in  Ni-based  superal¬ 
loys  [14],  and  a  quantitative  model  for  the  creep  rate  for 
microtwinning  was  proposed  [21].  The  energetics  associ¬ 
ated  with  the  reordering  process  was  explored  using  first 
principles  calculations  [22]  where  the  most  likely  diffusional 
paths  were  identified.  This  analysis  suggests  that  the  activa¬ 
tion  energy  barrier  for  reordering  should  be  very  similar  to 
that  for  Ni  self-diffusion  in  the  ordered  phase,  under  the 
simplifying  assumption  of  Ni3Al  being  representative  of 
the  Y  structure.  In  addition  to  controlling  the  rate  of  ther¬ 
mally  activated  microtwinning,  it  has  been  postulated  that 
reorder-mediated  shearing  may  also  be  rate  limiting  to  the 
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Fig.  1.  Schematic  illustration  of  twinning  in  the  Ll2  structure  of  the  y'  precipitate  phase,  (a)  Creation  of  a  true  twin  structure  by  passage  of  a/ 3(1  1 2)  super 
Shockley  partial  dislocations  on  adjacent  {111}  planes;  (b)  creation  of  a  pseudotwin  structure  of  a/ 6(1 1 2)  Shockley  partial  dislocations  on  adjacent  {111} 
planes. 
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formation  of  both  superlattice  intrinsic  and  extrinsic  stack¬ 
ing  faults  in  superalloys  at  elevated  temperatures  [22], 
Based  on  the  concept  of  the  generalized  stacking  fault 
energy  landscape,  phase  field  modeling  has  been  used  to 
parametrically  study  microstructural  effects  on  the  criti¬ 
cally  resolved  shear  stress,  which  was  also  coupled  with 
temperature  effects  on  the  atomic  reordering  kinetics 
related  to  dislocation  glide  [23].  The  aforementioned  rate 
limiting  process  helps  explain  the  dynamic  motion  of  twin¬ 
ning  partial  dislocations  and  microtwin  propagation  cap¬ 
tured  by  Legros  et  al.  [17]  during  elevated  temperature 
in  situ  TEM  straining  experiments  of  a  rafted  y/y'  superal¬ 
loy,  where  it  was  observed  that  the  twinning  partial  dislo¬ 
cations  moved  at  a  higher  velocity  in  the  y  phase  and  a 
lower  velocity  in  the  y'  phase. 

Recent  research  on  creep  deformation  mechanisms  in 
polycrystalline  Ni-based  turbine  disk  superalloy  Rene  104 
has  shown  that  microtwinning  is  a  dominant  deformation 
mechanism  during  creep  at  a  test  temperature  of  677  °C 
and  a  stress  level  of  690  MPa  [16],  Deformation  in  this  case 
was  carried  out  to  such  a  high  level  of  plastic  deformation 
(2.0%  strain)  that  it  was  difficult  to  capture  discrete  loca¬ 
tions  where  microtwins  may  have  initiated  or  terminated. 
Due  to  the  present  lack  of  experimental  evidence  on  micro¬ 
twin  initiation  and  evolution  there  has  yet  to  be  a  detailed 
description  of  how  the  microstructure  and  other  external 
variables  (stress  and  temperature)  dictate  the  microtwin¬ 
ning  process  from  initiation  to  propagation.  Therefore, 
the  motivation  for  this  study  was  to  perform  an  in-depth 
experimental,  theoretical,  and  computational  study  on  the 
evolution  of  the  microtwinning  deformation  mode  by 
focusing  on  the  salient  microstructural  features  that  pro¬ 
mote  microtwinning  over  other  deformation  mechanisms. 

2.  Materials  and  experimental  procedure 

The  material  used  in  this  study  is  a  powder  metallurgy 
product,  polycrystalline  Ni-based  turbine  disk  alloy 
(Rene  104)  that  was  supplied  in  the  form  of  a  scaled 
up  turbine  disk  forging  by  General  Electric  Aircraft 
Engines  (GEAE,  Evendale,  OH).  The  nominal  composi¬ 
tion  of  this  alloy,  which  has  Co,  Cr,  Al  and  Ti  as  the 
major  alloying  elements,  is  given  elsewhere  [24].  The  disk 
forging  was  subjected  to  a  proprietary  two  step  supersol- 
vus  heat  treatment,  which  refers  to  a  final  heat  treatment 
wherein  the  primary  y'  precipitates  that  are  present  dur¬ 
ing  the  powder  compaction,  extrusion,  and  forging  pro¬ 
cesses,  are  re-solutionized  back  into  the  matrix,  then 
precipitated  as  a  multimodal  y'  precipitate  size  distribu¬ 
tion  [25,26].  For  the  material  in  this  study  the  resulting 
microstructure  consisted  of  a  bimodal  size  distribution 
of  Y  precipitates,  with  the  coarser  being  the  secondary 
Y  and  finer  being  the  tertiary  y'  precipitates.  On  average 
the  secondary  and  tertiary  y'  precipitate  sizes  (and  vol¬ 
ume  fractions)  were  experimentally  measured  to  be 
186.6  (48.4%)  and  25.5  nm  (2.7%),  respectively,  using 
quantitative  image  analysis  on  material  that  was 


Fig.  2.  Microstructure  of  Ni-based  disk  superalloy  Rene  104  following  a 
supersolvus  heat  treatment.  SEM  image  shows  the  bimodal  y'  precipitate 
size  distribution  (average  secondary  and  tertiary  y'  precipitate  sizes  186 
and  25  nm,  respectively). 

extracted  from  the  rim  section  of  the  turbine  disk  forging 
(Fig.  2).  Specimens  for  creep  testing  and  TEM  deforma¬ 
tion  mechanism  identification  were  extracted  by  wire 
electric  discharge  machining  (EDM)  from  the  same  loca¬ 
tion  and  prepared  as  standard  round  tensile  creep 
specimens. 

In  order  to  track  the  evolution  of  the  microtwins  and 
to  determine  the  salient  microstructural  features  that  give 
rise  to  microtwinning  the  tests  were  interrupted  at  vary¬ 
ing  stages  of  creep  deformation.  Constant  load  tensile 
creep  experiments  were  preformed  at  677  °C  and 
724  MPa  and  interrupted  at  nominally  0.05%,  0.2%, 
and  1.0%  strain.  Specimens  were  cooled  under  load  to 
preserve  the  deformation  substructure.  Following  creep 
testing  thin  foils  for  TEM  substructure  characterization 
studies  were  extracted  from  the  gage  at  a  45°  angle  with 
respect  to  the  tensile  axis  in  order  to  maximize  the  prob¬ 
ability  of  imaging  deformation  activity  along  slip  planes 
that  experienced  maximum  shear  stress.  The  foils  were 
mechanically  thinned  using  progressively  finer  SiC  paper 
to  a  thickness  of  150  pm  then  slurry  drilled.  Final  thin¬ 
ning  was  performed  in  a  Struers  Tenupol  5  twin  jet  pol¬ 
ishing  unit  using  an  electrolyte  consisting  of  60% 
methanol,  35%  2-n-butoxyethanol,  and  5%  perchloric 
acid  at  a  temperature  of  45  °C  and  applied  voltage  of 
15  V.  A  detailed  diffraction  contrast  TEM  analysis  was 
performed  using  a  Philips  CM200  or  a  FEI  Tecnai 
TF20  microscope,  operated  at  200  keV.  High  resolution 
TEM  (HRTEM)  and  high  angle  annular  dark  held 
(HAADF)  imaging  were  performed  in  a  probe-corrected 
FEI  Titan  held  emission  gun  transmission  electron 
microscope  operated  at  300  keV. 
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3.  Results 

3.1.  Microtwinning  as  the  dominant  deformation  process  at 
larger  strain 

Following  creep  at  677  °C  to  0.98%  plastic  strain  over 
2725  h  microtwinning  was  found  to  be  the  dominant  mode 
of  deformation.  Fig.  3a  and  b  shows  a  bright  held  (BF) 
TEM  image  of  the  microtwins  (viewed  nearly  edge  on 
and  oriented  along  the  [Oil]  zone  axis)  and  the  corre¬ 
sponding  selected  area  diffraction  (SAD)  pattern,  respec¬ 
tively.  The  SAD  pattern  shows  fundamental  reflections 
from  the  y  matrix,  superlattice  reflections  from  the  y'  pre¬ 
cipitates,  and  the  [111]  twin  reflections,  confirming  that 
the  planar  defects  are  microtwins.  Further  support  is  pro¬ 
vided  by  the  HRTEM  images  shown  in  Fig.  4,  where  twins 
of  several  different  thicknesses  are  shown  which  can  be  as 
narrow  as  two  {111}  planes  (Fig.  4a)  or  as  thick  as  several 


Fig.  3.  (a)  BF  TEM  micrograph  of  microtwins  formed  during  creep 
deformation  to  0.98%  strain  at  677  °C  and  724  MPa.  (b)  SAD  pattern 
showing  the  fundamental  y  matrix,  y'  precipitate  superlattice,  and  twin 
reflections  (marked  by  arrows).  B  =  [0 1 1], 


Fig.  4.  HRTEM  images  of  microtwins  present  within  the  y  matrix  phase 
at  different  thicknesses. 


tens  of  atomic  planes  (Fig.  4b  and  c).  The  two  layer  fault 
structure  is  recognized  as  an  extrinsic  stacking  fault  in 
the  matrix  and  as  a  superlattice  extrinsic  stacking  fault  in 
the  precipitates  [22,27]. 

Although  the  a/3  ( I  12)  partial  is  recognized  to  be  the 
twinning  partial  in  the  f  phase  [19],  extensive  investiga¬ 
tions  using  conventional  diffraction  contrast  imaging  [14] 
and  HAADF  scanning  transmission  electron  microscopy 
(STEM)  imaging  [22]  has  shown  that  microtwins  propagate 
by  the  motion  of  identical  a/6  ( I  12)  Shockley  partial  dislo¬ 
cations  that  shear  through  both  the  y  matrix  and  y'  precip¬ 
itates  on  adjacent  {111}  glide  planes,  an  example  of  which 
is  shown  in  Fig.  5.  As  mentioned  previously  and  discussed 
in  detail  elsewhere  [22],  shearing  of  the  precipitates  is 
believed  to  be  thermally  activated  due  to  diffusion-medi¬ 
ated  reordering  of  the  high  energy  faults  that  would  be  cre¬ 
ated  in  the  wake  of  the  Shockley  partials.  We  now  focus  on 
the  processes  that  occur  during  the  earlier  stages  of  creep 
deformation  that  ultimately  lead  to  this  microtwinned 
substructure. 

3.2.  Initiation  of  microtwins 

In  order  to  identify  the  nucleation  sources  of  the 
microtwins  it  was  necessary  to  examine  the  deformation 
substructure  at  the  earliest  stages  of  deformation.  TEM 
examination  following  0.05%  plastic  strain  revealed  a  par¬ 
ticular  grain  of  interest,  which  allowed  a  unique  view  of  the 
deformation  activity,  as  viewed  down  the  (111)  glide  plane. 
A  low  magnification  BF  STEM  image  of  this  particular 
grain  is  shown  in  Fig.  6.  The  grain  was  bounded  by  typical 
high  angle,  serrated  grain  boundaries,  as  well  as  an  anneal¬ 
ing  twin  boundary,  which  is  observed  on  the  lower  edge  of 
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Fig.  5.  Representative  HAADF  STEM  image  of  microtwin  shearing 
through  both  the  y  matrix  and  y'  precipitate. 


Fig.  6.  BF  STEM  image  of  the  deformation  substructure  following  creep 
to  0.05%  strain  at  677  °C  and  724  MPa  showing  deformation  activity 
being  emitted  from  second  phase  particles  and  grain  boundaries. 


the  image.  A  number  of  carbide  and/or  boride  phases  are 
present,  both  along  grain  boundaries  and  within  the  grain 
interiors.  These  second  phase  particles  are  typically  found 
in  this  and  other  polycrystalline  superalloys  produced  by 
powder  metallurgical  processes  [28].  For  instance,  the 
chemical  composition  of  the  intragranular  particles  as 
determined  by  energy  dispersive  X-ray  spectroscopy  tended 
to  be  enriched  in  Ti,  Ta,  Co,  Cr,  Mo,  and  Nb,  indicating  a 
MC  type  carbide  [29].  Both  M23C6  carbides,  enriched  in 
Co,  Cr,  Mo,  C  and  B,  as  well  as  borides,  have  also  been 
observed  along  grain  boundaries.  Upon  closer  examination 
these  minor  carbide  and  boride  phases  appear  to  be  distinct 
origins  of  deformation  activity.  The  nascent  deformation 
activities  near  these  second  phases  are  in  the  form  of  both 


dislocations  and  stacking  faults  in  the  matrix.  Note  that  the 
stacking  faults  viewed  in  this  orientation,  with  the  fault 
plane  nearly  parallel  to  the  foil  plane,  do  not  show  the 
alternating  bright  and  dark  fringes  normally  expected  when 
characterizing  stacking  faults.  Instead,  the  stacking  faults 
appear  as  dark  grey  regions  and  are  present  within  the  y 
matrix. 

3.3.  Detailed  TEM  deformation  substructure  analysis 

Additional  diffraction  contrast  TEM  analysis  was 
required  to  confirm  that  dislocations  emitted  from  the 
grain  boundary  carbides/borides  or  intragranular  carbide 
sources  readily  dissociate  into  Shockley  partial  disloca¬ 
tions.  Fig.  7  presents  a  series  of  higher  magnification  BF 
TEM  images  (using  different  two  beam  diffraction  condi¬ 
tions)  of  the  defect  structures  near  one  of  the  intragranular 
carbides  shown  in  Fig.  6.  The  specific  dislocation  configu¬ 
ration  highlighted  in  the  schematic  illustration  of  Fig.  7f 
indicates  that  the  leading  partial  dislocation  has  “perco¬ 
lated”  between  the  secondary  f  precipitates,  leaving  them 
unsheared,  while  some  of  the  finer  scale  tertiary  f  precipi¬ 
tates  have  been  looped  by  the  leading  partial.  The  stacking 
fault  and  associated  Shockley  partial  dislocations  were 
determined  to  be  on  the  (111)  glide  plane  (which  is  close 
to  the  foil  normal  in  this  grain)  based  on  the  g  •  R  invisibil¬ 
ity  criterion  for  planar  faults  in  the  fee  crystal  structure, 
where  R  is  the  fault  vector.  The  stacking  fault  is  visible  in 
Fig.  7a  since  g  =  Ill  and  g  •  R  is  not  zero.  For  this  case 
the  leading  partial  dislocation  is  visible  but  the  trailing 
one  is  not.  Conversely,  using  g  =  11 1  (Fig.  7b)  the  trailing 
dislocation  is  visible  while  the  leading  one  is  not.  Addi¬ 
tional  g  vectors  were  chosen  such  that  g  ■  R  =  0  and  the 
contrast  due  to  the  partial  dislocations  was  not  obscured 
by  the  stacking  fault.  For  instance,  the  trailing  partial  dis¬ 
location  is  visible  for  g  =  220  (Fig.  7c)  and  g  =  022 
(Fig.  7d),  while  it  is  invisible  when  imaged  with  g  =  202 
(Fig.  7e).  The  leading  partial  dislocation  is  visible  when 
imaged  with  g  =  220  and  g  =  202,  while  it  is  invisible  with 
g  =  022.  Based  on  this  analysis,  the  leading  partial  disloca¬ 
tion  has  a  Burgers  vector  of  a/6  [2  1  1]  and  the  trailing  par¬ 
tial  dislocation  has  a  Burgers  vector  of  a/6[I  2  I],  produced 
by  the  dissociation  of  a  perfect  a/ 2[1  1  0]  matrix  disloca¬ 
tion.  For  a  clearer  illustration  of  this  dislocation  configura¬ 
tion  a  schematic  is  presented  in  Fig.  7f  with  regard  to  the 
(111)  plane  of  the  Thompson’s  tetrahedron.  In  addition 
to  determining  the  Burgers  vectors  of  the  dislocations,  it 
is  also  evident  that  a  single  slip  system  was  operative  in  this 
region,  since  all  the  dislocations  can  be  accounted  for  by 
contrast  from  these  perfect  and  Shockley  partial 
dislocations. 

With  increasing  creep  deformation  to  0.2%  strain  the 
same  stacking  faults  are  observed;  however,  they  are  pres¬ 
ent  at  a  much  higher  density.  The  local  plastic  zone  near  an 
intragranular  carbide  is  shown  at  a  higher  TEM  magnifica¬ 
tion  in  Fig.  8,  which  shows  several  a/ 2(1  10)  dislocations 
originating  from  an  intragranular  carbide  source.  In  this 
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Fig.  7.  A  series  ol'  BF  TEM  images  showing  the  diffraction  contrast  analysis  of  perfect  and  partial  dislocations  bounding  a  stacking  fault  in  the  y  matrix 
using  the  following  imaging  conditions:  (a)  g  =  Ill,  B  =  [2  1 1];  (b)  g  =  ill,  B  =  [1 1 0];  (c)  g  =  220,  B  =  [1  1  1];  (d)  g  =  022,  B  =  [1 1  1];  (e)  g  =  202, 
B  =  [1  1  1],  (fj  Schematic  illustration  of  the  dislocation  configurations  along  with  Thompson  tetrahedron  notation  of  the  Burgers  vector. 


case  dislocation  emissions  on  inclined  {111}  glide  planes 
are  observed.  The  dislocations  directly  adjacent  to  the  car¬ 
bide  are  comprised  predominantly  of  narrowly  dissociated 
a/2(110)  dislocations,  all  contained  within  the  y  matrix. 
There  are  also  a  few  instances  where  significantly  wider  dis¬ 
sociation  events  occurred,  creating  wide  stacking  faults 
(labeled  F  in  Fig.  8a).  Fig.  8b  shows  the  defect  structure 
a  few  hundred  nanometers  from  the  same  carbide,  exhibit¬ 
ing  primarily  stacking  faults  in  the  matrix,  with  Shockley 
partials  looping  the  secondary  y'  precipitates.  Thus,  these 


observations  suggest  that  while  a/ 2(1 10)  dislocations  are 
initially  generated  near  carbides,  a  progression  towards  a 
fully  dissociated  (decorrelated)  condition  occurs,  where 
the  Shockley  partials  comprising  the  c//2(  l  1 0)  dislocations 
are  widely  separated. 

It  is  important  to  note  that  these  stacking  fault  struc¬ 
tures  appear  to  be  present  at  a  similar  density  to  that 
observed  for  microtwins  formed  at  larger  strain,  which  fur¬ 
ther  confirms  the  importance  of  the  former  as  a  precursor 
to  microtwinning.  In  addition,  nearly  all  of  the  grains 
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Fig.  8.  BF  TEM  images  of  (a)  the  plastic  zone  surrounding  an 
intragranular  carbide,  showing  emission  of  dissociated  glide  a/2(110) 
dislocations,  with  a  few  dissociated  into  wider  stacking  faults  (marked  F) 
and  (b)  widely  dissociated  a/ 2(1 1 0)  dislocations  in  a  region  approximately 
1  pm  from  the  carbide. 


examined  via  TEM  possess  these  deformation  features. 
When  the  faults  are  either  inclined  or  nearly  edge  on  with 
respect  to  the  foil  normal,  it  is  observed  that  a  majority 
of  the  faults  do  not  shear  through  the  y'  precipitates.  As 
mentioned  in  the  previous  section,  this  indicates  that  the 
decorrelated  Shockley  partial  dislocations  bounding  the 
stacking  faults  do  not  shear  through  the  precipitates  them¬ 
selves,  since  this  would  leave  a  complex  stacking  fault. 
Instead,  the  «/6(112)  Shockley  partial  dislocation  loops 
that  are  left  around  the  precipitates  have  the  same  Burgers 
vector  as  the  leading  Shockley  partial  dislocation,  and  an 
intrinsic  stacking  fault  is  left  within  the  matrix. 

The  nature  of  the  stacking  faults  at  0.2%  creep  strain 
were  determined  through  complementary  BF/DF  TEM 
imaging  pairs  under  a  two  beam  diffraction  condition  with 
the  deviation  parameter  set  close  to  s  =  0  and  the  operative 
diffraction  vector  g  =  1 1 1  and  g  =  TIT.  As  seen  in  Fig.  9a 
and  b,  the  outer  fringes  in  the  BF  image  are  symmetrical 
and  the  outer  fringes  in  the  centered  DF  image  are  asym¬ 
metrical.  Several  of  the  stacking  faults  overlap  in  this 


Fig.  9.  (a)  BF  TEM  and  (b)  centered  DF  TEM  images  of  stacking  faults  in 
the  y  matrix  developed  during  creep  deformation  to  0.2%  strain  at  677  °C 
and  724  MPa. 


image,  so  this  analysis  is  focused  on  the  three  stacking 
faults  between  the  two  secondary  y'  precipitates  in  the  cen¬ 
ter  of  the  image  in  Fig.  9.  To  determine  whether  the  stack¬ 
ing  faults  are  intrinsic  or  extrinsic  in  nature  the  g  vector  is 
placed  in  the  center  of  the  stacking  fault  in  the  centered  DF 
image  (Fig.  9b).  Under  these  imaging  conditions  the  g  vec¬ 
tor  points  away  from  the  bright  outer  fringe,  which 
indicates  that  the  stacking  faults  are  intrinsic  based  on 
two-beam  dynamical  theory  [30],  Thus  the  stacking  fault 
structures  formed  in  the  matrix  are  similar  in  nature  to 
those  observed  at  smaller  strain  levels. 

4.  Discussion 

4.1.  Analysis  of  dislocation  decorrelation 

The  concept  of  decorrelation  of  matrix  dislocations  due 
to  interaction  with  precipitates  has  been  discussed  recently 
by  Raujol  et  al.  [31]  and  is  based  on  their  experimental 
observations  in  both  turbine  disk  and  turbine  blade  alloys. 
Moreover,  the  experimental  evidence  described  above  indi¬ 
cates  that  decorrelation  of  a/ 2(1 1 0)  matrix  dislocations  is  a 
necessary  precursor  to  microtwinning.  While  the  possible 
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mechanisms  for  twin  thickening  will  be  discussed  below,  in 
this  section  we  focus  on  developing  a  fundamental  under¬ 
standing  of,  and  a  predictive  model  for,  the  conditions 
for  dislocation  decorrelation. 

4.1.1.  Phase  field  simulation 

In  order  to  analyze  all  possible  dislocation  activities  for 
a  given  y  channel  width  as  a  function  of  the  orientation  and 
magnitude  of  the  in-plane  shear  stress  phase  field  simula¬ 
tions  of  dislocation-precipitate  interactions  were  carried 
out.  Incorporating  the  generalized  stacking  fault  (GSF) 
[32]  energy  of  the  octahedral  slip  plane,  the  phase  field 
method  [33,34]  naturally  describes  dislocation  dissociation 
and  the  associated  stacking  faults  and  partial  dislocations 
in  both  the  ordered  and  disordered  phases.  The  results 
(Fig.  10)  show  that  a  rich  variety  of  dislocation  activities 


can  be  observed,  including  Orowan  looping  and  anti-phase 
boundary  (APB)  shearing  by  a  full  dislocation,  Orowan 
looping  by  both  partial  dislocations,  Orowan  looping  by 
only  the  leading  partial,  and  stopping  of  both  partials  at 
the  entrance  to  the  y'  channel.  Periodic  boundary  condi¬ 
tions  create  an  array  of  spherical  secondary  f  particles 
with  a  radius  of  150  nm.  The  y  channel  width  between 
the  y'  particles  is  77.1  nm,  as  shown  in  the  schematic  draw¬ 
ing  (lower  right  corner  of  Fig.  10).  The  Burgers  vector 
direction  for  the  perfect  a/ 2(1 1  0)  dislocation  is  along  the 
vertical  direction  (set  as  0°),  thus  the  Burgers  vector  for 
the  leading  partial  AS  is  along  +30°  and  the  trailing  partial 
SB  is  along  -30°,  as  illustrated  schematically  in  the  center 
of  Fig.  10.  The  angle  6  between  the  in-plane  shear  stress 
direction  and  the  Burgers  vector  direction  of  the  perfect 
dislocation  AB  can  vary  from  —180°  to  +180°  (however, 
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Fig.  10.  Dislocation  activity  diagram  comparing  the  analytical  formulae  with  the  phase  field  simulation  results  (lines  are  from  the  analytical  formulae, 
symbols  are  calculated  from  phase  field  simulations).  (A)  Orowan  bowing  of  full  dislocation;  (B)  Orowan  bowing  of  both  partial  dislocations  (dislocation 
dissociation);  (C)  APB  shearing  of  y';  (D)  Orowan  bowing  of  the  leading  partial  dislocation  (decorrelation);  (E)  dislocations  arrested  in  front  of  the 
channel.  Note  that  the  leading  partial  Burgers  vector  is  AS,  while  that  for  the  trailing  partial  is  SB. 
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for  6  >  90°  or  0  <  —90°,  SB  should  be  defined  as  the  leading 
partial  and  the  dislocation  should  initially  be  placed  on  the 
right  side  of  the  y'  particles).  As  indicated  in  the  plot,  each 
dot  in  Fig.  10  corresponds  to  one  phase  field  simulation 
performed  under  that  specific  stress  magnitude  and  orien¬ 
tation  0.  Different  symbols  represent  different  dislocation 
activities,  which  are  explained  in  the  inset  table  in  Fig.  10. 


4.1.2.  Analytical  modeling 

Based  on  the  phase  field  simulation  results,  analytical 
expressions  have  been  developed  to  determine  the  condi¬ 
tion  for  a/ 2(1 1 0)  dislocation  dissociation  as  a  function  of 
the  magnitude  and  orientation  of  the  shear  stress  i,  the 
intrinsic  stacking  fault  (ISF)  energy  y1  ,  the  friction  force 
on  partial  dislocations  tf  due  to  solute  drag  or  cutting  of 
smaller  tertiary  y'  precipitates  [35],  and  the  influence  of 
microstructure  (y  channel  width  between  secondary  y'  par¬ 
ticles).  This  development  is  based  on  the  classical  Orowan 
bowing  mechanism  [36],  For  simplicity  the  dependence  of 
the  line  tension  on  the  screw/edge  character  of  the  disloca¬ 
tion  (i.e.  the  angle  between  the  dislocation  sense  vector  and 
the  Burgers  vector)  is  ignored,  thus  the  line  tension  is 
assumed  to  be  T  =  (ctpb2)/2,  where  a  is  a  constant  whose 
typical  value  is  1-2. 

Depending  on  the  orientation  of  the  Burgers  vector  of  the 
dislocation  with  respect  to  the  applied  stress,  the  leading 
Shockley  partial  could  have  a  higher  or  lower  resolved  shear 
stress  compared  with  the  trailing  one.  In  general,  decorrelat- 
ed  motion  of  the  partial  dislocations  can  take  place  only  if  the 
leading  partial  has  a  higher  resolved  shear  stress.  Exceptions 
may  exist  considering  the  reduction  in  the  effective  channel 
width  for  the  trailing  partial  to  pass  through  when  the  lead¬ 
ing  one  deposits  segments  in  the  channel,  which  repel  the 
trailing  one,  especially  for  narrow  y  channel  widths  and 
low  ISF  energy  cases,  which  will  be  elaborated  later. 

The  detailed  formulation  of  the  analytical  model  is  pre¬ 
sented  in  the  Appendix  A.  According  to  the  model  the  crit¬ 
ical  in-plane  stress  for  the  leading  partial  to  bow  through 
the  channel  t[rss  (which  corresponds  to  the  configuration 
shown  in  Fig.  11a)  is  expressed  as: 

j  2tf  a  fib 

Tcrss  cos(0i)  +  cos(0t)  +  2cos(0i) 
for  cos(0i)  >  cos(0t) 

j  2  /  a  f.ib  \ 

Tcrss  cos(0i)  +  cos(0t)  \  2  tfJ 

for  cos(0|)  ^  cos(0t)  (1) 


TcrssCOS(0t)  =  A  + 


a  ■  [.ib 
2-2  H 


yISF 

~b~ 


(2) 


T[rss  is  the  corresponding  critical  stress  for  the  trailing  par¬ 
tial  dislocation  to  pass  through  the  channel,  H  is  the  dis¬ 
tance  between  the  deposited  leading  particle  segments  in 
the  channel  and  the  trailing  dislocation,  which  is  a  function 
of  r*rss,  and  yISF  is  the  intrinsic  stacking  fault  energy. 

These  analytical  equations  of  the  critical  stresses  are  val¬ 
idated  against  the  phase  field  simulation  results.  Materials 
parameters  used  in  both  the  phase  field  simulations  and  the 
analytical  formulae  are  as  follows:  {111}  plane  shear  mod¬ 
ulus  for  both  the  y  and  y'  phases  is  p  =  67.2  GPa,  the  ISF 
energy  is  lOmJm-  ,  and  the  dislocation  line  tension  is 
T  =  2.45(/(/r/2).  Note  that  the  line  tension  chosen  for  these 
calculations  is  relatively  high  (i.e.  a  =  2.45,  which  is  higher 
than  the  typical  a  value  of  1-2  as  previously  mentioned) 
because  the  extra  gradient  energy  term  in  the  phase  field 
dislocation  model  at  the  length  scale  indicated  in  Fig.  10 
leads  to  a  relatively  high  line  tension  [37],  In  order  to  make 
a  quantitative  comparison  with  the  phase  field  simulation 
results  the  same  line  tension  is  used  in  the  analytical  model. 
In  addition,  the  numerical  grid  in  the  phase  field  simula¬ 
tions  effectively  creates  a  friction  force  (Peierls  stress)  cali¬ 
brated  as  Tf=100MPa  and  the  same  friction  force  is 
assumed  in  the  analytical  model  as  well.  Considering  the 
diffuse  nature  of  dislocation  cores  in  the  phase  field  simula¬ 
tions  at  such  a  coarse  grained  length  scale  the  channel 
width  for  the  leading  dislocation  is  77.1  nm  and  that  for 
the  trailing  one  is  65.4  nm.  The  same  channel  widths  are 
assumed  in  the  analytical  calculations. 

The  analytical  result  (solid  red  line  and  dashed  blue  line) 
is  superimposed  on  the  phase  field  simulation  results  in 
Fig.  10.  The  solid  red1  curve  is  the  critical  stress  for  the 
leading  dislocation  to  pass  through  the  channel  and  the 
dashed  blue  curve  is  the  critical  stress  for  the  trailing  dislo¬ 
cation  to  pass  through.  As  observed,  the  analytical  calcula¬ 
tions  agree  reasonably  well  with  the  phase  field 
simulations.  The  deviation  of  the  critical  stress  for  the  lead¬ 
ing  partial  at  high  0  angles  is  an  artifact  that  will  be  dis¬ 
cussed  in  the  following  section.  As  will  be  shown  below, 
implementation  of  the  analytical  formulae  allows  a  quick 
identification  of  the  conditions  (stress  magnitude,  orienta¬ 
tion,  intrinsic  stacking  fault  energy,  y  channel  width,  fric¬ 
tion  force,  etc.)  under  which  decorrelated  motion  of 
Shockley  partials  occurs  for  a  given  microstructure,  as  well 
as  helps  avoiding  issues  associated  with  coarse-graining  to 
the  length  scale  considered  in  this  paper  such  as  high  dislo¬ 
cation  line  tension  and  high  friction  force. 


where  0j  is  the  angle  between  the  applied  stress  and  the  Bur¬ 
gers  vector  of  the  leading  partial  and  0t  is  the  angle  between 
the  stress  and  the  Burgers  vector  of  the  trailing  partial.  The 
channel  width  is  2  and  tf  is  the  friction  force  caused  by  sol¬ 
ute  drag  or  the  cutting  of  tertiary  particles. 

The  force  balance  at  the  critical  moment  for  the  trailing 
partial  dislocation  to  bow  through  the  channel  (which  cor¬ 
responds  to  the  configuration  shown  in  Fig.  lib)  is: 


4.1.3.  Effect  of  applied  stress  on  dislocation  decorrelation 
From  the  applied  tensile  stress  in  the  experiment 
(724  MPa)  the  maximum  in-plane  shear  stress  is  calculated 
as  362  MPa,  which  is  indicated  in  Fig.  12  by  the  boundary 


1  For  interpretation  of  color  in  Figs.  1,  4,  7,  10-14,  the  reader  is  referred 

to  the  web  version  of  this  article. 
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Fig.  11.  Schematic  illustrations  depicting  (a)  the  critical  moment  for  the  leading  Shockley  partial  to  pass  through  the  channel  and  (b)  the  critical 
configuration  for  the  trailing  Shockley  partial  dislocation  to  pass  through  the  y  channel. 


m>a  o  +15 


Fig.  12.  Critical  stress  plot  (friction  stress  if  =  50  MPa,  channel  width  64  nm,  yIS  =10  mJ  m~2,  line  tension  T  =  (\xb2)/2.  The  solid  red  line  (marked  L)  is 
the  critical  stress  for  the  leading  partial  calculated  using  the  analytical  form  from  Eq.  (1);  the  solid  blue  line  (marked  T)  is  the  critical  stress  for  the  trailing 
partial  calculated  using  the  solution  of  Eq.  (2);  the  dotted  red  line  (marked  L')  is  the  critical  stress  for  the  leading  partial  calculated  using  a  simplified  form; 
the  dotted  blue  line  (marked  T'J  is  the  critical  stress  for  the  trailing  partial  calculated  using  a  simplified  form. 


of  the  shaded  area.  The  materials  parameters  used  in  the 
calculations  were:  elastic  constants  of  the  y  matrix, 
Cn  =  224.3  GPa,  C12  =  148.6  GPa,  C44  =  125.8  GPa  [38]; 
the  corresponding  shear  modulus  for  the  {111}  slip  plane 
u  =  67.2  GPa;  the  Poison’s  ratio  v  =  c'l.  =  0.4;  the  lat- 
tice  constant  of  the  y  matrix  phase  a  =  0.358  nm;  the  Bur¬ 
gers  vector  of  the  Shockley  partials  b  =  0.146  nm.  The 


friction  force  due  to  solute  drag  or  shearing  of  tertiary  par¬ 
ticles  is  estimated  to  be  50  MPa  [3 1  ].  Since  this  estimation 
was  made  for  room  temperature  it  is  expected  to  be  an 
upper  bound  of  the  friction  force  operative  at  the  creep 
temperature  of  677  °C.  Although  the  experimental  results 
indicate  that  the  intrinsic  stacking  fault  (ISF)  energy  yISF 
is  ~25  30  mJ  m~2  at  intermediate  temperatures  [39],  a 
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more  recent  quantitative  analysis  of  the  Shockley  partial 
dissociation  suggests  that  yISF  %  10  niJ  m“2  [31].  Thus 
y1  =  10  mJ  nT  2  was  chosen  for  the  calculations.  For  sim¬ 
plicity  the  line  tension  was  assumed  to  be  T  =  (jilr/2) 
regardless  of  its  character  (edge  or  screw).  Fig.  12  shows 
the  calculation  results  obtained  for  a  y  channel  width  of 
64  nm,  which  corresponds  well  with  the  average  channel 
width  measured  experimentally.  The  red  curve  (marked 
L)  is  the  critical  stress  for  the  leading  partial  to  pass 
through  the  channel,  while  the  blue  curve  (marked  T)  is 
the  critical  stress  for  the  trailing  partial  to  pass  through. 
It  can  be  seen  from  Fig.  12  that  when  the  resolved  shear 
stress  on  the  leading  partial  becomes  smaller  than  that  on 
the  trailing  one  (0  <  0)  the  two  curves  collapse,  indicating 
no  decorrelated  motion. 

For  a  shear  stress  level  of  362  MPa,  as  long  as  the  angle 
6  between  the  shear  stress  direction  and  the  full  dislocation 
Burgers  vector  direction  is  between  +70°  and  +33°,  decor- 
related  motion  of  partials  is  expected  because  the  leading 
partial  can  go  through  the  64  nm  channel  while  the  trailing 
partial  cannot.  For  angles  from  +33°  to  the  branching 
point  of  the  curve  (marked  X  on  the  plot,  around  +9°) 
both  dislocations  can  go  through  the  channel,  although 
the  leading  one  still  has  a  higher  resolved  stress,  which 
means  it  moves  faster  than  the  trailing  one.  In  this  case  dec¬ 
orrelated  motion  is  still  possible,  but  depends  on  the  differ¬ 
ence  between  the  critical  stress  for  the  leading  and  trailing 
partials.  For  0  between  +9°  to  —53°  the  partials  will  pass 
through  the  channel  together  as  the  trailing  one  pushes 
the  leading  one,  and  thus  no  decorrelated  motion  will  occur 
and  shearing  by  perfect  a/2(\  1 0)  dislocations  should  dom¬ 
inate  the  deformation  process.  For  9  between  —53°  and 
—90°  no  dislocation  is  able  to  pass  through  the  channel 
at  this  stress  level. 

The  decorrelation  diagram  with  respect  to  stress  orienta¬ 
tion  has  mirror  symmetry  relative  to  the  normal  of  the 
a/2(110)  full  Burgers  vector  (the  dislocation  will  move 
from  the  right  to  the  left  in  this  case,  because  of  reversal 
of  the  in-plane  shear  stress  orientation).  Hence,  there  will 
be  large  asymmetry  for  a  given  slip  vector  depending  upon 
the  sign  of  the  applied  shear  stress.  This  asymmetry  is 
brought  about  by  the  reversed  role  of  the  leading  and  trail¬ 
ing  Shockley  partials.  Moreover,  for  a  {111}  slip  plane 
there  are  three  co-planar  a/ 2(1 10)  dislocations,  AB,  BC, 
and  AC,  whose  Burgers  vectors  are  120°  apart.  By  consid¬ 
ering  the  dissociation  of  all  three  types  of  full  dislocations 
the  probability  of  observing  decorrelated  dislocation 
motion  for  a  random  orientated  362  MPa  in-plane  shear 
stress  is  about  78%  for  an  average  channel  width  of  64  nm. 

Douin  et  al.  [40]  also  proposed  a  set  of  simple  analytical 
forms  to  describe  the  dissociation  process  that  provides  the 
critical  stress  of  the  leading,  trailing,  and  full  dislocations. 
The  results  are  plotted  in  Fig.  12  (incorporating  the  addi¬ 
tional  solute  drag  friction  stress  if)  as  dashed  lines  (the  line 
for  leading  dislocation  is  marked  L',  the  line  for  trailing 
dislocation  is  marked  T',  and  the  line  for  full  dislocation 
is  marked  S,  within  the  —90°  <  9  <  90°  region).  Clear  dis¬ 


agreements  exist  between  the  two  calculation  results.  Com¬ 
pared  with  Eq.  ( 1 )  the  critical  stress  calculated  for  line  L' 
does  not  take  into  account  the  assistance  from  the  trailing 
dislocation  through  elastic  interaction.  Thus  Eq.  ( 1 )  gives  a 
lower  critical  stress  when  9  C  60°.  However,  when  9  >  60° 
(for  —90°  <  9  <  90°,  when  AS  is  still  defined  as  the  “leading 
partial”)  the  resolved  shear  stress  on  the  trailing  partial  SB 
becomes  negative,  which  means  that  SB  is  driven  to  move 
in  the  opposite  direction  with  respect  to  AS.  In  this  case 
the  scenario  shown  in  Fig.  11a  is  no  longer  accurate 
because  the  trailing  partial  tends  to  move  away  from  the 
leading  partial.  Under  this  condition  the  critical  stress  cal¬ 
culated  for  line  L',  which  is  established  for  large  separation 
distances  between  the  partials,  can  better  describe  the  crit¬ 
ical  stress  for  the  leading  partial.  As  for  the  trailing  dislo¬ 
cation,  the  critical  stress  calculated  for  line  T'  does  not 
take  into  account  the  reduction  in  effective  channel  width 
due  to  the  segments  laid  by  the  leading  partial,  which  in 
turn  underestimates  the  critical  stress  level. 

4.1.4.  Effect  of  microstructure  on  dislocation  dissociation  and 
microtwinning 

Previous  work  has  shown  that  for  deformation  condi¬ 
tions  that  produce  microtwins  in  the  finer  scale  precipitate 
microstructure  considered  in  this  work  a  different  deforma¬ 
tion  mode  appears  to  be  dominant  for  a  coarser  precipitate 
microstructure  developed  by  a  slower  cooling  rate  from  the 
supersolvus  temperature  [41].  For  this  latter  microstruc¬ 
ture,  which  simulates  that  in  the  “bore”  region  of  a  turbine 
disk,  deformation  occurs  by  intense  a/2(l  10)  dislocation 
activity  in  the  matrix,  with  occasional  shearing  of  the  pre¬ 
cipitates  by  superlattice  intrinsic  or  extrinsic  stacking 
faults.  Thus  for  this  coarser  microstructure  dislocation  dec¬ 
orrelation  does  not  seem  to  occur.  Rather,  matrix  perfect 
dislocations  seem  to  retain  their  ability  to  percolate 
through  the  y  channels.  The  analytical  model  developed 
above  provides  an  insight  into  this  microstructure  effect. 

In  Fig.  13  the  average  channel  width  is  assumed  to  be 
159  nm,  which  corresponds  to  the  microstructure  observed 
in  the  bore  region.  As  expected,  the  wider  channel  makes  it 
easier  for  dislocations  to  pass  through,  and  the  critical  stress 
level  decreases.  When  compared  with  Fig.  12,  it  can  also  be 
observed  that  the  intersection  of  the  critical  stress  for  the 
trailing  dislocation  to  pass  through  the  channel  with  the 
362  MPa  applied  shear  stress  contour  is  shifted  from  +33° 
to  +50°  as  the  channel  width  increases.  Thus,  decorrelation 
of  partials  is  less  likely  (in  terms  of  the  angular  range  of  the 
shear  stress).  Considering  all  three  co-planar  a/2(l  1 0)  dislo¬ 
cations,  the  probability  of  observing  decorrelated  disloca¬ 
tion  motion  for  a  randomly  oriented  in-plane  shear  stress 
of  362  MPa  is  51%  for  a  channel  width  of  1 59  nm,  compared 
with  78%  for  a  channel  width  of  64  nm,  considered  earlier. 
(Note  that  the  probabilities  considered  here  are  for  decorre¬ 
lation  of  any  of  the  three  co-planar  slip  systems  under  a  ran¬ 
domly  oriented  362  MPa  in-plane  shear  stress.)  Thus, 
depending  on  the  applied  stress,  there  will  be  a  certain  range 
of  channel  widths  that  maximize  the  probability  of 
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Fig.  13.  Critical  stress  plot  (friction  stress  if  =  50  MPa,  channel  width  159  nm,  yISF  =  10  mJ  m~2,  a  =  1,  line  tension  T  =  (|i62)/2.  The  solid  red  line  is  the 
critical  stress  for  the  leading  partial  calculated  using  the  analytical  form  from  Eq.  (1);  the  solid  blue  line  is  the  critical  stress  for  the  trailing  partial 
calculated  using  the  solution  of  Eq.  (2). 


decorrelated  dislocation  motion.  For  instance,  if  the  in-plane 
stress  is  200  MPa  instead  of  362  MPa,  then  for  a  channel 
width  of  64  nm  even  the  leading  partial  dislocation  cannot 
pass  though.  For  a  channel  width  of  159  nm  dissociation  is 
still  possible.  Using  these  analytical  formulae  we  may  iden¬ 
tify  the  optimal  stress  level  for  a  given  microstructure,  or 
specify  target  microstructures  for  a  given  applied  stress,  that 
will  favor  matrix  dislocations  or  microtwinning  as  the  dom¬ 
inant  deformation  mode. 

4.2.  Transition  from  stacking  faults  to  microtwins 

A  comprehensive  review  of  deformation  twinning  has 
been  provided  by  Christian  and  Mahajan  [19].  Possibilities 
for  the  creation  of  microtwins  from  intrinsic  matrix  stack¬ 
ing  fault  structures  have  been  identified.  The  pole  mecha¬ 
nism  has  been  proposed  to  describe  twinning  in  fee 
crystals  [42].  In  this  model  it  is  envisioned  that  a  twin  is 
formed  from  a  perfect  a/ 2(1  10)  dislocation  that  acts  as 
an  anchor  point  where  the  twinning  dislocation  (a/ 6(1 1 2) 
Shockley  partial)  sweeps  out  on  {111}  glide  planes  and 
rotates  about  the  pole  to  produce  the  twin.  Alternatively, 
in  the  model  proposed  by  Mahajan  and  Chin  [43]  two 
co-planar  a/ 2(1 1 0)  dislocations  with  different  Burgers  vec¬ 
tors  are  envisioned  to  react  and  form  three  identical  Shock- 
ley  partial  dislocations  on  adjacent  {111)  glide  planes, 
thereby  acting  as  the  nucleus  of  a  three  layer  twin.  From 
the  limited  experimental  evidence  of  microtwinning  in 
polycrystalline  f  strengthened  superalloys  investigated 
recently  [14-16]  it  appears  that  all  of  the  a/ 2(1 1 0)  disloca¬ 
tions  that  are  observed  to  dissociate  into  Shockley  partial 
dislocations  have  the  same  Burgers  vector.  Accordingly, 
these  observations  suggest  that  the  microtwinning  mecha¬ 
nism  in  these  alloys  differs  from  that  proposed  by  Mahajan 
and  Chin.  Furthermore,  in  the  double  cross-slip  model  pro¬ 
posed  by  Pirouz  a  single  screw  oriented  a/ 2(1 10)  disloca¬ 
tion  dissociates  into  two  Shockley  partial  dislocations 
[44],  Similarly  to  a  Frank-Read  source,  the  large  shear 
stress  imposed  on  the  leading  Shockley  partial  dislocation 
will  bow  out  and  create  an  intrinsic  stacking  fault.  As  the 


leading  Shockley  partial  recombines  with  the  static,  trailing 
partial  the  a/ 2(1 1 0)  dislocation  that  is  once  again  created 
can  cross-slip  (along  its  screw  orientation)  to  an  adjacent 
{111}  plane,  then  redissociate  into  Shockley  partial  dislo¬ 
cations  in  the  same  manner  as  the  dissociation  on  the  initial 
{111}  plane.  If  this  operation  were  to  be  repeated  on  mul¬ 
tiple  adjacent  {111}  planes  then  Shockley  partial  disloca¬ 
tions  possessing  the  same  Burgers  vector  would  be 
operative,  and  a  microtwin  could  potentially  develop. 
Based  on  the  present  work  it  is  impossible  to  definitively 
validate  this  mechanism,  however,  it  is  noted  that  this 
mechanism  is  consistent  with  the  operation  of  single 
u/2(l  1 0)  slip  systems  as  the  source  of  the  twins. 

The  present  observation,  that  extended  intrinsic  faults 
form  in  the  y  matrix  preceding  the  onset  of  microtwinning, 
brings  to  bear  another  possibility  for  microtwin  formation. 
It  has  been  shown  that  these  intrinsic  faults  form  from  the 
dissociation  of  u/2(110)  dislocations  that  originate  near 
intragranular  carbides  and  intergranular  carbides/borides. 
Under  appropriate  shear  stress  conditions  the  matrix  dislo¬ 
cations  undergo  decorrelation,  with  the  leading  partial  dis¬ 
location  percolating  between  the  secondary  f  particles  and 
laying  out  large  areas  of  intrinsic  stacking  faults.  With 
these  intrinsic  stacking  faults  established,  the  subsequent 
interaction  of  another  matrix  dislocation  will  create  a  sce¬ 
nario  such  as  that  illustrated  in  Fig.  14.  Here  it  is  envi¬ 
sioned  that  the  same  a/ 2(1 10)  dislocation  type  cross-slips 
onto  the  intrinsic  stacking  fault  plane,  then  dissociates. 
As  the  leading  partial  of  this  second  dislocation  drives  for¬ 
ward  it  creates  an  extrinsic  stacking  fault  in  the  matrix.  As 
this  partial  reaches  locations  where  the  first  leading  partial 
has  been  arrested  at  f  particles,  the  pair  of  identical  Shock- 
ley  partials  are  then  able  to  traverse  f  via  reorder-mediated 
shearing  [22].  This  process  will  result  in  the  removal  of  the 
Shockley  partial  loops  surrounding  the  f  particles,  as 
observed  for  instance  in  Fig.  7,  and  the  creation  of  SESF 
in  the  sheared  f  particles  following  reordering.  Repetition 
of  this  process  of  matrix  dislocation  cross-slip,  extension, 
and  reorder-mediated  shearing  would  enable  thickening 
of  the  extrinsic  stacking  faults/SESF  into  multilayer 
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Fig.  14.  Possible  scenario  for  transition  from  dissociated  a/ 2(1 10)  type  dislocations  and  decorrelated  leading  (Lj)  and  trailing  (Tj  Shockley  partial 
bounding  matrix  stacking  faults  into  microtwins  through  a  sequence  of  cross-slip  and  dissociation  steps  of  subsequent  a/ 2(1  1 0)  dislocations.  After  /  is  (a) 
looped  by  the  leading  Shockley  partial,  the  key  steps  are  (b  and  c)  the  approach  and  cross-slip  of  a  second  a/ 2(1  10)  dislocation,  (d)  dissociation  of  the 
second  a/ 2(1 10)  dislocation  (D2)  on  the  {111}  plane  adjacent  to  the  stacking  fault,  (e)  cooperative  shearing  of  y'  by  two  Shockley  partials  by  reorder- 
mediated  motion,  creating  superlattice  extrinsic  stacking  in  y'  and  an  extrinsic  stacking  fault  in  the  matrix,  (f)  approach  of  a  third  a/ 2(1 10)  dislocation 
(D3).  Repetition  of  this  process  would  enable  thickening  into  microtwins. 


microtwins.  This  scenario  is  also  consistent  with  the  obser¬ 
vation  of  microtwins  in  the  absence  of  significant  second¬ 
ary  slip  system  activity  [45],  Clearly,  the  trailing  Shockley 
partials  that  do  not  participate  in  the  forward  shearing  pro¬ 
cess  would  accumulate  near  the  origin  of  the  twins.  Their 
mutual  stress  fields  could  eventually  inhibit  further  thicken¬ 
ing  of  the  microtwins  by  this  cross-slip  process.  This  may 
provide  an  explanation  for  the  predominance  of  thin 
microtwins,  with  few  as  thick  as  that  shown  in  Fig.  4c. 

5.  Conclusions 

Through  a  detailed  TEM  characterization  study  on  the 
early  stages  of  creep  deformation  in  a  y'  strengthened  Ni- 


based  superalloy  it  was  possible  to  identify  discrete  loca¬ 
tions  where  deformation  microtwins  may  have  nucleated. 
Microtwinning  is  an  important  creep  deformation  mode 
at  low  stress  and  intermediate  temperatures  (650-750  °C). 
It  was  found  that  carbides  and  borides,  which  are  either  sit¬ 
uated  at  the  grain  boundaries  or  within  the  grain  interiors, 
acted  as  sources  for  a/ 2(1 1 0)  matrix  dislocations.  At  small 
strain  levels  matrix  dislocations  are  frequently  observed  to 
dissociate  (i.e.  decorrelate)  into  Shockley  partial  disloca¬ 
tions,  creating  broad  intrinsic  stacking  faults  in  the  y 
matrix.  Under  similar  creep  conditions  but  larger  strain 
levels  microtwinning  is  observed  to  dominate  the  deforma¬ 
tion.  Thus  matrix  dislocation  decorrelation  appears  to  be  a 
key  precursor  to  microtwinning. 
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The  decorrelation  process  is  driven  by  the  relative  forces 
acting  on  the  two  Shockley  partial  dislocations  of  a  perfect 
a/ 2(1 1 0)  -type  matrix  dislocation  in  the  narrow  confines  of 
the  y  matrix  channels,  since  the  y  precipitates  cannot  be 
sheared  by  the  leading  Shockley  partials.  A  detailed  analy¬ 
sis  of  these  forces  has  been  performed  using  computer  sim¬ 
ulations  based  on  a  phase  field  dislocation  model,  upon 
which  an  analytical  model  was  also  developed.  In  the  anal¬ 
ysis  a  consistent  picture  emerges  demonstrating  the  possi¬ 
bility  of  a  rich  variety  of  behaviors  depending  on  several 
factors,  including  orientation  of  the  applied  shear  stress  rel¬ 
ative  to  the  dislocation  Burgers  vectors,  the  y  channel 
width,  and  the  matrix  stacking  fault  energy.  The  models 
predict  that  decorrelation  can  occur  at  low  stress  levels  rel¬ 
ative  to  other  possible  processes  (such  as  Orowan  bowing 
around  particles)  when  the  applied  shear  stress  is  suitably 
oriented  so  as  to  create  a  larger  force  on  the  leading  Shock- 
ley  partial  dislocation.  This  result  provides  a  natural  expla¬ 
nation  for  the  observation  of  matrix  stacking  faults 
dominating  the  substructure  of  many  grains  at  small  strain 
levels.  Several  possible  scenarios  for  the  subsequent  thick¬ 
ening  from  intrinsic  stacking  faults  into  microtwins  have 
also  been  discussed,  including  repeated  cross-slip  and  dis¬ 
sociation  of  matrix  dislocations  onto  pre-existing  stacking 
faults. 
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stress  on  the  leading  particle  /rss  for  passage  through  the 
y  channel  between  two  y'  particles  is  given  by  the  equation: 

y  +  Tcrssc°s(0i)- A  =  °^  (A1) 

where  1  is  the  channel  width.  Assuming  that  the  trailing 
segment  close  to  the  channel  has  a  local  curvature  of  t\ 
(see  arrow  in  Fig.  11a),  the  force  balance  on  the  trailing 
one  is: 


^crss  COS(0t)-Tf- 


a  •  /(b 
2  rt 


£_ 

b 


(A2) 


The  local  curvature  of  the  trailing  partial  at  this  critical 
moment  can  be  estimated  if  we  assume  that  the  resolved 
shear  stress  is  much  greater  than  the  friction  force,  i.e. 


T 

crss 

a  •  f. ib 
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(A3) 


Thus 


2rt  =  X 


cos(0i) 

cos(0t) 


(A4) 


when  cos(0i)  >  cos(0t).  When  cos(0j)  <  cos(0t)  the  re¬ 
solved  shear  stress  on  the  leading  partial  is  lower  than  or 
equal  to  that  on  the  trailing  one  and  hence  the  trailing  par¬ 
tial  will  push  the  leading  one.  In  this  case  the  two  disloca¬ 
tions  will  have  similar  curvatures  at  the  critical  moment 
(indicated  in  Fig.  1  la  by  the  dotted  line  for  the  trailing  dis¬ 
location),  i.e. 


2  rt  «  X 


(A5) 


When  t\  is  known  the  critical  in-plane  stress  for  the  lead¬ 
ing  partial  to  bow  through  the  channel  is: 

j  2tf  a  jib 

Tcrss  cos(0j)  +  cos(0t)  +  Xcos(6\) 

forcos(0i)  >  cos(0t) 

1  2  ( \ 

Tcrss  cos(0j)  +  cos(0t)  V  ^  / 

for  cos(0i)  <  cos(0t)  (A6) 


Appendix  A.  Analytical  model  of  dislocation  dissociation 

During  the  critical  moment  when  the  leading  partial  is 
attempting  to  bow  through  a  y  channel  between  two  y'  par¬ 
ticles,  as  indicated  in  Fig.  11a,  the  resolved  shear  stress  on 
the  leading  partial  is  tcos(0[)  and  that  on  the  trailing  par¬ 
tial  is  t cos(0t),  where  0j  is  the  angle  between  the  applied 
stress  and  the  Burgers  vector  of  the  leading  partial  and  0t 
is  the  angle  between  the  stress  and  the  Burgers  vector  of 
the  trailing  partial,  as  shown  in  Fig.  12.  If  the  angle 
between  the  stress  and  the  Burgers  vector  of  the  full  dislo¬ 
cation  is  0,  then  0t  =  30°  +  0  and  0]  =  30°  —  0. 

The  leading  partial  experiences  a  repulsion  force/1  from 
the  elastic  interaction  with  the  trailing  partial  and  a  drag 
from  the  intrinsic  stacking  fault  between  them.  According 
to  the  Orowan  equation  [36]  the  critical  resolved  shear 


When  the  trailing  dislocation  attempts  to  pass  through 
the  channel,  no  matter  whether  the  leading  one  has  already 
passed  through  or  is  just  moving  together  with  and  right  in 
front  of  the  trailing  one  when  cos(0i)  <  cos(0t),  the  leading 
partial  must  have  already  laid  two  segments  at  the  y/y' 
interface  inside  the  channel  as  illustrated  in  Fig.  1  lb.  Thus 
the  effective  width  of  the  channel  is  reduced  for  the  trailing 
one  due  to  the  repulsive  force  of  the  deposited  segments  of 
the  leading  partial.  The  distance  H  between  the  deposited 
leading  particle  segments  in  the  channel  and  the  trailing 
dislocation  can  be  calculated  using  the  equation: 


cos  (Qt)b  +  yISF 


f.ib2 

4tl7/(1  —  v) 


(A7) 


which  represents  the  force  balance  between  the  resolved 
shear  stress,  the  attractive  force  from  the  intrinsic  stacking 
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fault  and  the  repulsive  force  from  the  leading  partial  seg¬ 
ment.  T‘rss  is  the  critical  resolved  shear  stress  for  the  trailing 
partial  to  pass  through  and  v  is  the  Poisson’s  ratio. 

Considering  that  the  effective  channel  width  for  the  trail¬ 
ing  partial  is  reduced  to  X  —  2 H,  the  force  balance  at  the 
critical  moment  for  the  trailing  partial  to  bow  through 
the  channel  is: 

•  ub  vISF 

<rss  COS(gt)  =  Tf  +  x  _  2H  ~  ~Y  (A8) 

As  this  is  a  quadratic  function  of  t'rss  (note  that  H  is  a 
function  of  T*rss)  an  analytical  solution  can  be  easily 
obtained.  The  other  condition  required  for  decorrelated 
motion  is  that  r*rss  7s  <TSS.  If  t',.ss  <  <TSS  the  trailing  partial 
will  push  the  leading  one  passing  through  the  channel  more 
or  less  simultaneously  and  thus  there  is  no  decorrelated 
motion. 
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